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Abstract 

In this work we use the lattice regularization method to study the behavior of the six point renor¬ 
malized coupling constant dehned at zero momentum for the three-dimensional theory in the 
intermediate and strong coupling domain. The result is in good agreement with the corresponding 
study in the Ising limit. 
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1 Introduction 


We study in this paper the behavior of the six point renormalized coupling constant dehned at 
zero momentum for the {(p‘^)d=‘i theory for a large variation of the bare coupling constant in the 
symmetric phase. For consistent we also show the behavior of the four point renormalized coupling 
constant. In such regime the perturbation theory doesn’t work and the nonperturbative approach 
is mandatory. As an analytic alternative approach in studying this theory in the strong coupling 
regime it is used the strong coupling expansion series [1, 2, 3]. The basic idea of this method 
consists in factoring out the kinematical parts of the Lagrange from the path integral and the 
result of evaluating the remaining non-Gaussian in close form is a formal expansion of the vacuum 
functional as a series in inverse power of coupling constant. It is possible from this expansion to 
extract a set of simple diagrammatic, which can be used to compute the n-point Green’s functions 
of the theory. For the {(p'^)d= 2 , theory it can be shown from the analytic study in the strong coupling 
limit the renormalized coupling constant approaches to the asymptotic limit. The same thing can 
be said for the six point renormalized coupling constant. From the point of view of numerical 
study a relative little work has been done in order to simulate the general theory in the domain 
of intermediate to strong coupling constant. For this reason we were motivated to prepare this 
paper. This work is a continuation of previous paper [4] where the six point renormalized coupling 
constant behavior of the (v?"^)d =2 theory was analyzed. The other important reason is that 
theory admits a non-trivial continuum limit [3, 5, 6] contrary to {(p^)d=i- Much more effort has 
been done in the application of renormalization group and high temperature expansion to the 
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theory {(p^)d =3 in the Ising limit [7]. The techniques that have been used to determine the proper 
n{n < 8)-point Green functions (IPI) include high-temperature lattice expansion, Monte Carlo 
methods and e expansion. For the case of d = 3 a complete list of references is given in [8]. Bender 
and Boettcher [9] did the analytic study for renormalized sextic coupling. This study takes into 
account the strong coupling calculation on hypercubic lattice in d-dimensions, which were invented 
to obtain the continuum limit. In two dimensions Sokolov and Orlov used renormalization group 
expansion and PadCBorel-Leroy resumation technique to get qq [10]. In the literature there are 
a few works available of the sextic coupling constant using the approach of lattice Monte Carlo. 
The hrst work was performed by Wheater [11] and later in the Ising limit by Tsypin [8]. So it 
seems appropriate the nonperturbative study using the lattice Monte Carlo technique [12] to get 
not only the asymptotic value of the higher order renormalized coupling constant but also the 
behavior on the intermediate and strong regime. Here as usual we consider Ti = c = 1. 


2 Higher order coupling constants 


We consider in the continuum the 0^ theory in d-dimensions Euclidean space in the presence of a 
source J, which the bare action is given by: 

^[0, J]= J (fx 

We introduce the vacuum persistent functional Z[J]\ 




( 1 ) 


Z[J]=J T>[0] exp(-^[0, J]) 


( 2 ) 
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and from this we define the generating functional W[J] for the connected Green functions, by 
writing 

W[J] = In Z[J]. (3) 


The vacuum expectation value of the held 0c is given by 


0c =< 0 > = 


5J{x 


-W[J] 


j=o 


(4) 


and the connected n-point Greens functions Gn{xi... Xn) is obtained from the generating func¬ 
tional (3), 


Gjii^Xi . . . Xyi) 
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5J{xi) ' ' ' 5J{Xn) 


InZ(J) 


(5) 


j=o 


The effective action r[0c] is dehned by a functional Legendre transform of hL[J]: 


r|0j = J d‘‘x Mx) J(x) - w\.j]. 


( 6 ) 


It is well known that the effective action is the generating functional of one-particle-irreducible(lPI) 
vertices, in particular the functional r[0] has a Taylor expansion in powers of 0 at 0 = 0; 


r[0c] f d'^Xi. . .d'^XnT^'"\xi, . . . ,Xn)4>c{Xl) . . .4>c{Xn)- 

„ n J 


(7) 


Here ... , Xn) is the proper n-point Green functions (IPI). Now we consider a source J which 

is constant and uniform in Euclidean space-time. This implies that 0 is a constant independent 
of space-time. We dehne the Fourier transform of P*^"'^(a;i,..., Xn) by 


T'-'‘\xu---,Xx)- j '■■■ j 


(2ir) 


dkn 

W) 




2=1 


( 8 ) 
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Thus, for constant 4>{x) = 0 we have 


r[0c] = E ^ 0, 0,, 0 ) 0 " (27r)''5(0) ^ (27r)"5(0) [/(0) (9) 

which is the defining equation for the effective potential U (0). Thus, the Taylor coefficients of the 
effective potential are the IPI vertices r*^"')(0, 0, 0,..., 0) evaluated at zero external momentum. 
The renormalized fourier transform proper n-point Green functions is obtained as fol¬ 

lows: the wave-function renormalization is obtained from the Fourier transform connected Green 


function of two-points from: 


1 ^ dG2{p^) 
dp^ 


and the renormalized mass is defined by 


ml = ZG^^ip^) 


In general the renormalized r^"')(0) proper n-point Green’s functions are given by 


rW(o) = ^ 


and it follows that the renormalized effective potential can be written as. 


G. = E 


rW(o) 


where (Jy^ = Z 0. From here we see that L0 is the generating function of one particle irreducible 


renormalized Green’s function at zero external momentum on all legs. Since we are doing the study 
of the theory in the symmetric phase, we have = 0. The particular interest to us is the 
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renormalized coupling constant r(‘^^(0) and the renormalized sextic coupling constant (r(®^(0)), 
which can be expressed in terms of Fourier transform of connected Green functions as follows: 


rW ( o ) = -zHg^'{p^))‘g4p^) 


(14) 


and 

ffl(O) = -Z=(G2 -‘(p"))« (g,(p") - 10GJ(p")G2-‘(p"))^. (15) 

The quantities that will be extracted from lattice Monte Carlo simulation are the dimensionless 
renormalized zero momentum scattering amplitudes g 2 n dehned by, 


“ 77^2r^-r^d+d('2r^)! ‘ 


(16) 


3 Simulation results 

A discrete version of action related to Monte Carlo simulation can be written as 

s[(P,J]= + + + , (17) 

^ ^ X X 

where is a vector of length a in the positive p-direction. It’s convenient to work with dimen¬ 
sionless quantities when making numerical simulation. Thus we rescale in d-dimensions the held, 
the mass and the coupling constant according to 0 = ^ and g = gQa^~‘^ to make 

the lattice action independent of the lattice constant a. So we have the action as follows: 

A = \ + e/i) - 4>{x)f + ^ ^ (18) 

^ X,II ^ X X 
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As usually we impose periodic boundary condition on fields: 


0(n + L^) = (j){n) for all n , (19) 

We use the standard Metropolis algorithm combined with the Wolff single cluster method [13] 
which is used to avoid the trapping into meta stable states due to the underlying Ising dynamics. 
We use the cluster algorithm using the embedded dynamics for 0^ theory, according to the action, 
[14] 

*S*Ising -EE |0(a; + (l)(j){x) \ s{x + fi)s{x) , (20) 

X fi 

where s{x) = sign(0(x)). Statistical errors are evaluated taking into account the autocorrelation 
time in the statistical sample generated by the Monte Carlo simulation. For simulation it is used 
32^ lattice. We choose correlation length = l/'nr^ as 1 << .^^ << -h to minimize finite size effects. 
For numerical simulation purpose it is appropriated to work with the equivalence Langrange: 

5[0] = <f){x)(j){x + ef,) + J2 + E - l)^ (21) 

X^fl X X 

where 

^2 ^ 1 ^^ _ 2 ^ 9 = ^- ( 22 ) 

The most part of our effort it is to find the values of A and k for different values of the bare 
coupling constants g for a fixed value of the renormalized mass, which is found within few percent. 
As we vary the coupling constant to a greater valne we have to put more negative the bare 
mass Considering the translational invariance of the correlations functions, one can choose to 
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approximate the momentum derivation in eq. (10) by variation of G 2 {p^) across one lattice spacing 
and in one direction in order to calculate the renormalized mass 


ml = 


(LV 

r<0(o)2>-<|0(p)|2>i 

UttJ 

<|0(p)|2> 


(23) 


where here 0 is the Fourier transform of the held and p = (^,0) is the smallest available non-zero 
momentum. 

To hgure out g 4 and values from simulation we used lattice version of eqs.(14) and (15), 

G4(/) G2\p^) I < 0(0)4 > -3 < 0(0)2 >2 


^4(4!) = 




p2=o 


< 4>(oy >" C 


(24) 


j»(6!) = 


10G|(p2)Gj‘‘(p2) _ G8(p2)GJ»(p' 






2d 

r 


p2=o 


= lOgi 


2 < 0(0)® > -15 < 0(0)4 >< 0(0)2 > +30 < 0(0)2 >3 


< 0(0)2 >3 


(25) 


The expressions above are prohibited to use in the weak coupling constant regime due to the large 
statistical errors. However we notice that eqs. (24) and (25) are efficient in the intermediate and 
strong coupling regime where the statistical errors are reasonable [15]. We also notice that the 
statistical errors increase with value {L/^r)- The values in this simulation we hnd that for the 
renormalized coupling constant and sextic coupling constant are consistent with the prediction of 
analytic and numerical methods in the Ising limit. Our prediction for g^ = 2.03 ± 0.072 in the 
Ising limit is in good agreement with value obtained by Tsypin [8], gg = 2.05 ± 0.15. In hgure 
hg.(l) we present the result of simulations for a large range of the bare coupling constant for 
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ge and g^. In order to be consistent we also inclnde in fignre fig.(2) the resnlt of simulations of 
64^ lattice in two dimensions which was taken from previous work [4], We observe both g 4 ^ and 
ge, for two and three dimensions, approach to asymptotic constant value. This behavior is in 
according with numerical and analytical predictions. From both figures we also can see that the 
variation of both quantities are relative small considering the large variation on g. This means 
the behavior of effective potential in the strong coupling regime does not change significatively 
for large variation of the coupling constant. For the computation of renormalized higher orders 
fl' 2 n {n > 3) it is difficult due to large statistical errors even in the regime of strong coupling . To 
avoid this problem it is useful to look at methods, which calculate the connected Green’s functions 
directly. One of this method was discussed by Drumond et al [16, 17]. They write down for an 
action with a source term, the Langevin equation describing the stochastic evolution of the field 
in the theory. Differentiating both sides of the Langevin equation successively with respect to the 
source field gives a set of slave equations, which describe the stochastic evolution of the estimators 
of Grenn’s function of higher orders [18]. The interesting study would be the ( 0^)3 theory in the 
broken phase for finite coupling constant [19]. 

4 Conclusions 

In this paper we studied on the lattice the behavior of the renormalized sextic coupling for the 
theory at intermediate and strong coupling constant domain. We performed a large vari¬ 
ation, from intermediate to strong coupling, of the bare coupling constant and notice that the 



2.2 

2 

1.8 

fl'e ^-6 

1.4 

1.2 

1 

0.8 
9a 

0.6 

0.4 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 g 1 

Figure 1: The behavior of and with the bare coupling constant for d = 3 


proper six-point Green functions (IPI) ge evaluated at zero external momentum has a regular 
behavior. This quantity has a definite asymptotic value in the limit of 5 ^ ^ cx). Its qualitative 
behavior is similar to the two dimensions theory. We remark that in the expansion of effective 
potential in powers of the field the (j)^ term is not negligible. Finally our results for g 4 and g^ are 
in good agreement with the values available in the literature. 
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Figure 2: The behavior of and with the bare coupling constant for d=2 
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